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■ We propose a fast algorithm for computing the economic capi- 

. tal, Value at Risk and Greeks in the Gaussian factor model. The 

algorithm proposed here is much faster than brute force Monte Carlo 
simulations or Fourier transform based methods 1. While the algo- 
& . rithm of Hull- White [Hj is comparably fast, it assumes that all the 

loans in the portfolio have equal notionals and recovery rates. This is 
a very restrictive assumption which is unrealistic for many portfolios 
encountered in practice. Our algorithm makes no assumptions about 
[ the homogeneity of the portfolio. Additionally, it is easier to imple- 

ment than the algorithm of Hull- White. We use the implicit function 
theorem to derive analytic expressions for the Greeks. 
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1 The Gaussian Factor Model 



Let us consider a portfolio of N loans. Let the notional of loan % be equal to 
TVj. Then the zth loan represents fraction 

of the notional of the whole portfolio. This means that if loan i defaults and 
the entire notional of the loan is lost the portfolio loses fraction or 100/ £ % 
of its value. In practice when a loan i defaults a fraction of its notional 
will be recovered by the creditors. Thus the actual loss given default (LGD) 
of loan i is 

LGDi — fi(l — ri) (2) 

fraction or 

LGDi — 100/i(l - ri)% (3) 

of the notional of the entire portfolio. 

We now describe the Gaussian m-factor model of portfolio losses from 
default. The model requires a number of input parameters. For each loan i we 
are give a probability Pi of its default. Also for each i and each k — 1, . . . , m 
we are given a number w^k such that J2T=i w ik < 1- The number Wi^ is 
the loading factor of the loan i with respect to factor k. Let (pi, ... , (fi m and 
4>\i = 1, . . . , N be independent standard normal random variables. Let 
be the cdf of the standard normal distribution. In our model loan % defaults 
if 



W i,k4>k + 

fc=l 



\ 



k=i 



This indeed happens with probability p^. The factors (pi, ... , (fi m are usually 
interpreted as the state of the global economy, the state of the regional econ- 
omy, the state of a particular industry and so on. Thus they are the factors 
that affect the default behavior of all or at least a large group of loans in the 
portfolio. The factors 1 , . . . , <j) N are interpreted as the idiosyncratic risks of 
the loans in the portfolio. 
Let Ii be defined by 

1% I {loan i defaulted} (5) 

We define the random loss caused by the default of loan i as 

U = /■(! - r % )U, (6) 
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where is the recovery rate of loan i. The total loss of the portfolio is 

L = J> (7) 

i 

An important property of the Gaussian factor model is that the Lj's are 
not independent of each other. Their mutual dependence is induced by the 
dependence of each Li on the common factors (pi, ... , <p m . Historical data 
supports the conclusion that losses due to defaults on different loans are 
correlated with each other. Historical data can also be used to calibrate 
the loadings w^.the Lj's are not independent of each other. Their mutual 
dependence is induced by the dependence of each Lj on the common factors 
0i, ... , 4>m- Historical data supports the conclusion that losses due to defaults 
on different loans are correlated with each other. Historical data can also be 
used to calibrate the loadings w^. 



2 Conditional Portfolio Loss L 

When the values of the factors (f) m are fixed, the probability of the 

default of loan i becomes 

p * = $| ^(pO-^^w I (8) 



k 



The random losses Li become conditionally independent Bernoulli vari- 
ables with the mean given by 

E cmd (L i ) = fi(l-ri)p i (9) 

and the variance given by 

VAR cond (Li) = ff (1 - n) V (1 - p l ) (10) 

By the Central Limit Theorem [2] the conditional distribution of the port- 
folio loss L, given the values of the factors 4>u ■ ■ ■ , (fim, can be approximated 
by the normal distribution with the mean 

E con d{L) = ^2E cond (Li) (11) 
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and the variance 

VAR cond {L) = J2 VAR cond (Li) (12) 

i 

We define the CDF of the unconditional portfolio loss L by 

F(x) = Prob(L <x), (13) 

for all real numbers x. Since the conditional loss distribution is approximately 
normal, the CDF can be approximated by 

N 

F(x) w / <S>(E cond , VAR cond , x) J] p(0O#i • • • d(f> m , (14) 

where &(E cond , VAR cond , x) is the CDF of the normal distribution with mean 
E cond and variance VAR cond , while p(x) is the density of the standard normal 
distribution. 



3 Numerical Example 

In this section we apply the proposed algorithm to the single factor Gaussian 
model of a portfolio with N = 125 names. We choose a 125 name portfolio 
because it is the size of the standard DJCDX.NA.IG portfolio. We choose 
a single factor model because it is the one most frequently used in practice. 
The parameters of the portfolio are 



fi - 

Pi = 0.015 + 



1 

N 

0.05(i 



0.5 



N — 1 
O.l(z-l) 
N — 1 



0.1(2-1) , , 

= 0-5- J ^_ 1 ) , (15) 

where % — 1, . . . , N. 

In Figure ^ we compare the CDF computed using 10 6 Monte Carlo sam- 
ples with the CDF computed using formula ()14|1 . 1 The agreement between 

1 The author has the code implementing the algorithm described here in MATLAB, 
VBA for Excel and C. 
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the two is good. We add that the quality of approximation will improve even 
further for a bigger portfolio and many bank portfolios have much more than 
125 names. 



Figure 1: CDF of the Unconditional Portfolio Loss 
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4 Computation of VaR and Economic Capi- 
tal 

We define the VaR (value at risk) of a given portfolio as the level of loss 
(expressed as fraction of portfolio notional), such that the probability of the 
portfolio loss being less or equal to VaR is equal to a predefined confidence 
level q. Typically, q is chosen to be between 0.99 and 1. Thus we have 

Prob(L < VaR) = q (16) 
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From the results of the previous section it follows that an accurate approxi- 
mation to VaR can be found by solving the equation 

N 

J $(E cond , VAR cond , VaR) Y[ ■■■d(f) m = q (17) 

i=l 

This equation can be solved, for example, by bisection. To find the solution 
to (jT7jl with accuracy of 1 basis point (0.01%) we would need to evaluate the 
left hand side of ()17|) no more than 14 times. 

Example. We want to calculate the VaR of the portfolio (Q3J) with the 
confidence level q = 99.75%. We solve the equation fli7| ) and round the 
solution to the nearest basis point to arrive at VaR = 16.36%. We now run 
Monte Carlo simulation with 5, 000, 000 samples and compute the probability 
that the portfolio loss is less or equal to 16.36%. Rounded to the nearest basis 
point it turns out to be 99.75%. Thus the VaR of the portfolio is indeed 
16.36%. 

If desired, the convergence can be sped-up by using Newton's method. 
After the VaR is calculated we can calculate the economic capital by sub- 
tracting the average portfolio loss from VaR. 

We now compare the algorithm proposed here with other proposed al- 
ternatives. The FFT based methods pQ require the computation of a large 
number of Fourier transforms. To determine the VaR with an error of less 
than 1 basis point (0.01%) it is necessary to compute approximately 10,000 
Fourier transforms. Each Fourier transform is as expensive to evaluate as the 
left hand side of (fTTj) . Thus our algorithm is significantly faster than the FFT 
based methods. It is well known that the FFT methods are much faster than 
the direct Monte Carlo simulation. Thus our algorithm is much faster than 
the Monte Carlo approach. Finally, the recursive approach of Hull- White |3] 
is comparable in speed to the algorithm proposed here. However, it assumes 
that all the loans in the portfolio have equal notionals and recovery rates. 
This is a very restrictive assumption which is unrealistic for many portfolios 
encountered in practice. Our algorithm makes no assumptions about ho- 
mogeneity of the portfolio. Additionally, it is easier to implement than the 
algorithm of Hull- White. 
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5 Computing the Greeks 



Since VaR satisfies equation (JTTj) we can use the implicit function theorem 
[H] to find its partial derivatives with respect to the parameters of the model. 
These partial derivatives are traditionally called the Greeks in finance. We 
arrive at the following expressions 

d I ^{E cond ,VAR cond ,VaR)Uf =1 p{<Pi)d(pi . . .d<p m 

-VaR = - * — — „ /jw, T— ( 18 ) 



dN * I w^iHE cond , VAR cond , VaR) Uti 



d_ VaR _ I iiHEcond, VAR cond , VaR) Uli p(«Wi ■ ■ ■ #m 
d Pi a I W^R®( E cond, VAR cond ,VaR)U?=iP(<l>i)d<f>i...d<f) m 



d I d£~7®( Econd i VAR cond , VaR) Uf =1 p{4> i )d(j) 1 . . . 



VaR — — ^ ^ - __ _ (20) 



dw ^ I w^t®{Econdi VAR cond , VaR) Uti p(<k)#i ■ • • ■ 



d _ I j^jEcond, VAR cond , VaR) Y\li KWi • 

and 
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These expressions can be easily evaluated by numerical integration using 
Hermite-Gauss quadrature. 



6 Conclusions 

We proposed an algorithm for computing the VaR and the economic capital 
of a loan portfolio in the Gaussian factor model. The proposed method was 
tested on a portfolio of 125 names and gave high accuracy results. The 
accuracy will be even higher for portfolios with more names. Many of the 
bank portfolios are much larger than 125 names. The proposed algorithm 
is much faster than the FFT based methods and the brute force Monte 
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Carlo approach. The speed of the Hull- White algorithm |3] is comparable 
to that of the the algorithm proposed here, but the Hull- White algorithm 
requires that all the loans in the portfolio have equal notionals and recovery 
rates. This is a very restrictive assumption which is unrealistic for many 
portfolios encountered in practice. Our algorithm makes no assumptions 
about the homogeneity of the portfolio. Also, it is easier to implement than 
the algorithm of Hull- White. 

Additionally, we obtained analytical expressions for the Greeks using the 
implicit function theorem. 

We also comment that the algorithm can be extended trivially to the case 
of stochastic recovery rates and recovery rates correlated with the state of 
the factor variables. 

Some of the ideas used in this paper were previously explored in jl] and 
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